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SUMMARY 


In  this  report  the  stress  distribution  in  loop  and  bolt  joints  made  of 
compvjsite  materials  is  studied  using  analytical  and  experimental  methods. 
Also,  the  stress  concentration  in  ^^rthotropic  composite  sheets  with  a 
circular  hole  is  studied. 

The  numerical  technique  to  evaluate  the  stresses  and  displacements  is 
described.  Simple  design  formulas  to  compute  the  stresses  at  critical 
regions  are  developed.  A  procedure  of  structural  design  is  delineated  and 
applied  to  specific  cases. 

Experimental  results  obtained  by  testing  composite  joints  are  reported. 
Stress  distribution  tests  were  performed  using  the  photoelastic  technique. 
Failure  location  and  ultimate  loads  for  several  composite  joints  are  shown. 

This  report  contains  the  analysis  of  the  stress  concentration  that  appears 
around  a  circular  hole  in  an  orthotropic  cor>iposite  sheet  under  simple 
tension  and  simple  shear  loadings. 
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INTRODUCTION 


In  this  investij^nt  on,  methods  were  developed  for  calculating  and  design¬ 
ing  joints  of  comp»site  materials  that  have  directional  properties.  The 
introduction  of  loads  into  structures,  especially  into  thin  sheets,  his¬ 
torically  has  presented  problems,  even  in  isotropic  materials.  At  the 
beginning  of  this  century,  such  difficulties  were  overcome  for  the  case 
of  isotropic  materials.  Until  now,  however,  some  simplifying  assumptions 
for  joints  of  isot'opic  materials  have  been,  in  general,  tacitly  accepted. 


In  this  work,  for  ■■xample,  we  found  that  the  assumption  of  contact 
stresses  in  the  form  of  a  cosine  function  in  bolt  joints  yields  unreal¬ 
istic  results,  even  in  isotropic  materials.  For  the  present  investigation, 
the  contact  stresses  were  not  assumed  but  rather  were  an  outcome  contained 
in  the  method  of  solution. 


First,  a  stiffness  matrix  was  developed  by  which  composites  having  not 
only  one  but  also  n-directions  of  fiber  layers  could  be  taken  into 
account.  On  the  basis  of  this  development,  a  computer  program  was  estab¬ 
lished  through  which  the  stresses  and  displacements  in  a  joint  of  arbi¬ 
trary  shape  can  be  found.  To  illustrate,  three  sets  of  bolt  joints 
fabricated  from  glass  fiber  or  boron  fiber  composites  and  an  isotropic 
material  (e.g.,  steel)  were  calculated  and  studied  in  detail  for  several 
bolt  d iamoter -to- loop  thickness  relations.  Eighteen  bolt  joints  were 
manufactured,  tested,  and  compared  with  the  computer  results.  Finally, 
very  simple  relations  for  the  maximum  stresses  were  established. 


A  similar  approach  was  taken  to  the  calculation  of  stresses  in  cutouts  in 
composite  materials.  It  was  found  that,  in  the  contours  of  a  circular 
hole  cut  into  a  unidirectional  fiber  composite,  the  stress  concentrations 
are  about  twice  that  encountered  in  isotropic  materials  for  certain  types 
of  loading.  For  this  reason,  a  hole  in  a  fibrous  composite  must  be 
designed  so  that  tie  fibers  do  not  terminate  in  the  hole  but  rather  lead 
around  the  hole  in  an  uninterrupted,  organic  fashion.  If  this  is  done, 
then  cutouts  in  fibrous  composite  materials  will  be  highly  efficient  com¬ 
pared  to  their  counterparts  in  isotropic  materials. 


Stress  distribu 
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DISCUSSION 


METHOD  FOR  THE  STRESS  ANALYSIS  OF 
ANISOTROPIES  JOINTS  OF  ARBITRARY  FORM 

The  problem  is  to  find  che  stress  dis¬ 
tribution  at  a  joint  composed  of  an 
anisotropic  material  and  of  variable 
thickness  (Figure  1). 

In  general,  the  boundary  conditions 
are  mixed;  in  other  words,  displace¬ 
ments  or  stresses  are  given  along  the 
boundaries. 


To  find  the  solution,  we  divide  the 
joint  in  a  certain  finite  number  of 
triangles  and  assume  that  the  dis¬ 
placements  are  linear  functions  at 
each  one  of  these  small  triangles. 
Compatibility  and  equilibrium  equa¬ 
tions  for  each  node  of  this  irregular 
net  lead  us  to  a  system  of  algebraic 
equations  that  permit  us  to  obtain 
the  displacements  at  each  node  and, 
with  this,  to  determine  the  stress. 


Figure  1.  Subdivision  of  a 
Joint  Into  Finite 
Elements. 


Stiffness  Matrix  of  a  Triangular  Element 


For  a  triangular  element  (Figure  2),  we  assume  the  displacements  func- 

H  nn  ♦  »6  »8 

t  ion . 


U^  =  +  Sgl  +  SgTl 

Uy  =  +  ag?  +  agTl  (1) 

where  a.  are  constants. 

1 

The  displacements  of  the  nodes  1,  2,  and  3  are  given  by 


u  =  a, 

XI  ^ 
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Figure  2. 

Triangular  Element. 
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where  §  .  is  the  local  coordinate  system,  with  origin  at  point  1. 
In  matrix  form,  the  equations  (2)  are 
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or 


(M  =  [Al  {a] 

From  this  equation, 

{a}  =  (A)-'  [6] 

To  find  the  inverse  of  the  matrix  [A]  ,  we  write 

[A]  [a]  =  [I]  ,  [a]  =  [A]'" 

and  partition  in  the  following  manner: 
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Thus,  we  obtain  the  following; 
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Then,  inverting  th.e  matrix  [A*]  given  in  equation  (7),  we  obtain 
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Thus,  the  inverse  of  the  matrix  [A]  is 


[A] 


-1 


[A*J‘^  [0] 

(0]  [A*]‘^ 


where  [A*]  ^  is  given  by  equation  (9). 

From  equations  (1),  it  is  possible  to  obtain  the  strains 

Bu 

X 

X  a§  ^ 

9u 

3u  9u 
_ X  _ y  _ 

aTi  ^6 


xy 


(9) 

(10) 
(U) 


(12) 


(13) 


or,  in  matrix  form, 


0  1  0  0  0  0 

0  0  0  0  0  1 

0  0  10  10 


Symbolically,  we  write  equation  (14)  in  the  form 

[t]  =  [D]  •  {a}  (i: 

Substituting  {a}  given  by  equation  (5)  into  the  last  equation,  we 
obtain 

{e}  -  [D]  •  [A]‘M6}  (U 

On  the  other  hand,  the  relationship  between  stresses  and  strains  for  an 
anisotropic  material  and  the  two-dimensional  problem  is 
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The  strain  energy  is  given  by 


W  =  ^  r  It  a  +  €  a  +  y  t  idV 
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where  t  is  the  thickness  of  the  joint  and  {e)  is  the  transposed 
vector  of  (s}  .  But 

»  /(Dl  .  (a])''  -  {.f  (D)^ 

-  flAl"'  •  tM)''  iDl^  -  {6)^  /tAl"')^  ID)^  (21) 

where  we  used  equations  (15) and  (5).  With  equations  (19)  and  (21),  the 
strain  energy  is  expressed  by 

W  =  J  J  J  t  [D]^  [C]  [D]  [A]'^  {6}  dx  dy  (22) 


Now  we 


can  assume  that  the  forces  ^  >  Pg'ji  »  and 


P3T1  sre  applied  on  the  triangle  vertex,  as  shown  in  Figure  3. 

."’all 


1  P. 


Figure  3.  Generalised  Nodal  Forces. 


We  can  define  a  stiffness  coefficient  k*.  as  the  force  that  is  necessary 

1 1 

to  apply  at  i  direction  when  the  displacement  at  j  direction  is  the 
unit  and  all  other  displacements  are  zero.  The  directions  i  and  j  can 


be  any  one  of  Pj ^  , 
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XI  X3  yi 
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If,  for  example, 


X2 


1  and 


=  u  =  0  ,  the  stiffness  coefficients  k.  are  the 
ya  ya  12 


forces  that  we  need  to  apply  at  directions  of  Pi| 
keep  this  deformation  of  the  triangle. 


»  Pa  ( 


1-3^  ,  to 


To  compute  the  stiffness  coefficients,  we  will  use  the  second  Castigliano 
theorem,  which  permits  us  to  find  the  generalized  fo'ce  by  differentiation 
of  the  total  strain  energy  with  respect  to  the  corresponding  displacement. 
Thus , 


P.  = 
1 


bt. 


(23) 


Then,  by  differentiation  of  equation  (22)  with  respect  to  [6]  ,  we  have 

=  J  /  ^24) 

and  the  stiffness  matrix  (k]  yields 

Assuming  that  the  thickness  of  the  triangle  is  constant,  the  stiffness 
matrix  becomes 

[K]  =  ((Al''p  [D]'^  [C]  [D]  [A]’^  .  t  n  (26) 

with  n  equal  to  the  area  of  the  triangle. 

Equation  (24)  can  be  written  as 

{P}  =  IK]  {6]  (27) 

which  is  referred  to  the  axis  x,y  or,  in  other  words,  §  ,T1  .  However, 
for  certain  special  mixed  boundary  conditions,  it  would  be  useful  to 
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refer  ef|untion  (27)  to  the  rotated  axis.  The  rotations  are  Clj^  ,  Qtg  i 
and  tta  at  the  nodes  1,  2,  and  3,  respectively,  as  shown  in  Figure  4, 


Figure  4.  Rotation  of  the  Nodal  Directions. 
The  rotated  forces  vector  {P'}  is  related  to  {p}  as 


r 

Pix 

COS 

Ol 

0 

0 

-sin  tti 

0 

0 

r  '\ 

pi 

^ix 

Pax 

0 

cos  Og 

0 

0 

-sin  Ob 

0 

pi 

*^ax 
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pi 

Pax 

Ply 
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Ol 

0 
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ly 

Pay 

0 

sin  ttg 

0 

0 

cos  Og 

0 

Pay 

Pay 

k  ^  J 

0 

0 

sin  Oq 

0 

0 

cos 

^a 

P* 
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or,  in  condensed  form. 


An.ilogously ,  wt-  wrili' 


0  =  IT)  {S*} 


(30) 


By  substituting  equations  (29)  and  (30)  into  equation  (27),  we  obtain 


[P']  =  IK']  {6') 


(31) 


where  [K']  is  the  rotated  stiffness  matrix  given  by 


[K'l  =  (T]‘^  [K]  [T] 


(32) 


The  total  stiffness  matrix  of  tlie  whole  structure  can  be  obtained  by  sum¬ 
ming  in  adequate  order  the  elements  of  the  matrices  (K']  of  each  ele¬ 
ment.  The  law  of  formation  for  the  elements  of  the  total  stiffness 
matrix  [K*]  is 


N 


n=  I 


(33) 


where  k*  is  an  element  of  [K*]  and  k'.?  is  the  element  of  the 

i-J 

matrix  [K' ]  for  the  triangle  n  after  the  conversion  of  indices  i,j 
from  the  indices  1  through  6  to  the  actual  indices  of  the  whole  structure. 
The  summation  is  extended  to  the  N  triangles. 

If  the  numbering  of  the  nodal  directions  is  started  by  the  directions  in 
which  the  displacement  is  known,  the  matrix  equilibrium  equation  for  the 
whole  structure 


[K*]  {6*}  =  [P*3 


can  be  partitioned  as 


r 

"  * 

[Xu 

*  " 

K21 

• 

[k&] 

iin) 

►  * 

(p?) 

1 

(34) 


(35) 
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wlu'l 

re  {  61"  ] 

i  s 

tlie 

vec 

tor 

of 

i  s 

the 

vec 

tor 

of 

{prl 

i  s 

tlie 

vec 

tor 

of 

[P2] 

is 

the 

vec 

tor 

of 

The 

second  row 

of 

equ; 

ition  (35) 

[K^a] 

{ ^2 

■] 

the  known  nodal  displacements 
tlie  unknown  nodal  displacements 
the  unknown  nodal  forces 
tlie  known  nodal  forces 

can  be  written  as 

=  {pf}  -  [Kij  {6n 


(36) 


This  is  a  system  of  linear  equations  that,  when  solved,  enables  us  to 
know  the  displacements  [62}  .  Once  that  (62]  is  known,  the  first  row 
of  equation  (35)  gives  the  nodal  forces  {Pi^}  • 

Wlien  all  the  nodal  displacements  are  known,  the  stresses  can  be  determined 
by  substituting  the  displacements  {6}  ,  corresponding  to  each  triangle, 
in  equation  (19). 


It  must  be  noted  that  the  stresses  become  constant  at  each  triangle. 

Thus,  to  have  a  good  approximation,  the  triangular  net  must  be  more  dense 
in  those  places  where  the  stress  gradient  is  greater,  as  occurs  in  areas 
of  stress  concentrations. 


The  stress  diagrams  can  be  smoothed  out  adequately  by  taking  the  average 
of  stresses  between  two  adjacent  triangles  and  assigning  this  average  to 
the  midpoint  of  the  common  side. 


Computation  of  the  Matrix  [Cj  Elements 

In  the  most  general  case  we  can  assume  that 
the  material  contains,  at  an  arbitrary 
point,  fibers  in  n-directions.  For  the 
following  development,  we  take  only  three  y^ 
different  directions,  a,  b,  and  c,  but 
the  geneialization  to  any  number  of  dir¬ 
ections  is  immediate. 

We  refer  the  material  to  the  axes  x  and 
y  (Figure  5).  If  we  impose  a  deforma¬ 
tion  [e]  ,  the  materials  with  fibers  at 
directions  a,  b,  and  c  each  have  the 
stresses 


Figure  5.  Multilayered 
Composite . 


{^a)  =  [CgKe) 


[aj  =  [Cj(e]  (37) 
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where  the  vectors  {<^3}  1  '’i‘'’trices  [(’3]  , 

[Cjj]  ,  and  [C^]  are  referring  to  axes  x  and  y. 

Considering  a  certain  total  aiea  A  composed  of  the  areas  A^  ,  ^  > 

and  Aj.  corresponding  to  the  three  materials,  that  is  to  say, 

A  =  +  A^  +  A^  (38) 

the  forces  produced  by  the  imposed  deformation  are 

{Pal  =  A^  faj  =  A3  IC3]  [e] 

{Pbl  =  {^bl  =  \  [Cbl  {e} 

{Pel  =  {^cl  =  \  {Cj  {e}  (39) 

Then  the  total  force  is 

{p1  =  {Pfli  +  {Pbl  +  {Pel 


or,  by  using  equations  (39), 

{p1  =  ha  {^al  +  \  {^cl  )  {^1  (^0) 

On  the  other  hand,  it  is  possible  to  assume  that  [p]  is  produced  by  a 
stress  {a}  applied  on  the  total  area  A  , 

{P}  =  A  {a}  (41) 

Comparing  equations  (40)  and  (41),  we  find 

{c^l  =  (y  {^al  +  T  +  X 

Thus  we  arrive  at  a  matrix  [C]  given  by 

[C]  =  ng  [C3]  +  [C^]  +  n^,  [Cj  (43) 
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A  a  A|j 

wlR-re  ng  =  —  ,  "b  =  T  ’  "a  +  "b  +  "c  =  ^ 

If  the  matrix  ma  orial  (resin)  is  the  same  for  the  three  materials,  which 
is  tlie  most  coirni  jn  case,  the  numbers  ,  and  n^,  are  propor¬ 

tional  to  the  nu  iber  of  fibers  (or  layers)  at  each  direction  a  ,  b  , 
or  c . 

Referring  the  elastic  constants  for  the  material  a  to  the  axis  x  , 

7  ””  ^ 

Xg  (Figure  5) ,  we  have 


because  x^  ,  are  principal  axes  of  orthotropy  for  the  material  £  , 

Equation  (A5)  can  be  written  in  briefer  form,  as  follows: 

(a;)  =  [c;]  {€;)  (46) 

Now  we  want  to  transfer  equation  (46)  from  axis  Xg,ya  to  axis  x,y  . 

The  stresses  are  transformed  by  the  following  equation: 

=  (T13]  [aj  (47) 


where 

^a 

2 

■"a 

fTial  = 

2 

""a 

-2'^a"’a 

^a^a 

f2  2 

^3  * 

(48) 

and 

la  =  cos  cpa 

™a 

=  sin  cPa 

(49) 

The  strains  are 

transformed  by 

i^a) 

=  iT^al 

1 

(50) 
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with 


2 

"’a 

(-am^ 

iTaal 

= 

2 

"’a 

^a 

-^a"’a 

,2  2 
^a  -  "la 

(51) 


Substituting  equations  (47)  and  (50)  into  equation  (46)  and  then  multi¬ 
plying  both  members  by  (Tia  ]  ^  gives  us 


=  [Tial'"  (C^l  [Tg]  [s) 


(52) 


Thus,  by  comparison  of  equations  (52)  and  the  first  equation  (37), 
we  find 

[Cal  =  [Tia]*^  [C^]  [Tsai  (53) 

In  the  same  way,  it  is  possib’e  to  find  matrices  [Cj^]  and  [C^]  ,  using 
the  transformation  matrix  [T^bl  >  .  [Tab^  >  iTgcl  » 

Thus,  the  following  steps  are  necessary  to  find  matrix  [C]  of  a  triangle 

1.  Establish  angle  orientations  Cpa  »  <Pb  >  ^nd  of  the  fibers. 

2.  Establish  the  fiber  proportions  at  each  direction  n^  ,  nj^  , 
and  n^,  . 

3.  On  the  base  of  the  constants  for  each  layer,  obtain  [Ca]  , 

[Cjj]  ,  and  [Cj,]  by  using  equation  (53). 

4.  By  applying  equation  (43),  compute  [C]  . 

The  elastic  constants  for  each  layer  can  be  obtained  with  the  formulas 
for  unidirectional  composites  based  on  the  fiber  volumetric  content  and 
the  elastic  properties  of  the  fiber  and  matrix  material. 


Computer  Program 

With  the  elements  developed  in  this  section,  a  computer  program  was 
devised  which  enabled  us  to  establish  the  stresses  and  displacements  for 
a  joint.  The  shape  of  the  joint  can  be  arbitrary  as  well  as  the  system 
of  applied  loads.  In  other  words,  the  joint  can  be  solved  under  arbi¬ 
trary  mixed  boundary  conditions.  This  program  also  covered  a  very 
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specific  boundary  condition  that  appears  when  the  displacement  component 
at  the  normal  direction  of  a  cui-  fc?  boundary  is  restricted  and  also  when 
nonfriction  is  prescribed. 

The  thickness  of  the  joint  can  be  variable.  The  analysis  can  be  performed 
for  any  anisotropic  material  with  variable  elastic  properties. 


STRESS  ANALYSIS  OF  BOLT  JOINTS 


Basic  Assumptions 

The  stress  analysis  performed  during  this  contract  was  concerned  primarily 
with  finding  the  stresses  in  a  bolt  joint  under  tension  (Figure  6). 


Figure  6,  Typical  Joint. 


Between  the  bolt  and  the  joint  there  is  a  bushing  of  bronze  or  steel 
which  prevents  the  composite  material  from  coming  into  direct  contact 
with  the  steel  bolt. 

Critical  points  in  a  stress  analysis  of  this  type  of  joint  are  the  bound¬ 
ary  conditions  to  be  taken  along  the  hole  contour  in  contact  with  the 
bushing.  Many  investigators  adopted  the  cosine  law  for  the  pressure  dis¬ 
tribution  on  the  bolt.®  However,  the  stress  distribution  obtained  on  the 
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basis  of  tliis  assumption  does  not  correlate  with  the  experimental  results, 
especially  in  the  prediction  of  the  point  of  initial  failure.  In  fact, 
if  the  hypothesis  of  the  cosine  law  pressure  distribution  is  adopted,  the 
critical  point  is  A  in  I’igure  6  while  the  actual  point  is  indicated  at  a 
position  such  as  C.  A  more  realistic  stress  distribution  is  obtained  by 
assuming  that  the  radial  displacement  on  the  boundary  ACB  is  zero,  and 
that  there  is  no  friction  between  bolt  and  joint.  The  pressure  distri¬ 
bution  on  the  hole  boundary  obtained  from  this  hypothesis  shows  a  consid¬ 
erable  difference  with  respect  to  the  cosine  distribution.  Thus,  the 
boundary  conditions  adopted  in  the  following  analysis  are  those  depicted 
in  Figure  7. 


1.  On  the  external  boundary: 


2. 


3. 


nt 


0 


where  n,t  indicates  the 
normal  and  tangential  di¬ 
rections,  respectively. 

On  the  hole  boundary  ACB: 

u  =  0 

r 


with  r  and  9  ,  polar 
coordinates,  as  indicated 
in  yigure  7. 

On  the  hole  boundary  ADB: 


a 


r 


T 

re 


0 


These  boundary  conditions  correspond  to  a  tension  loading  on  the  joint. 


Stress  Distributions 


The  stress  distribution  of  joints  was  established  by  using  the  computer 
program  developed  on  the  basis  of  the  elements  given  in  the  preceding 
section.  A  parametric  study  was  conducted  to  find  the  differences 
between  the  stress  distribution  obtained  in  an  isotropic  joint  and  those 
in  fiber  reinforced  joints,  both  with  the  same  geometric  characteristics. 


16 


A  joint  was  selected  vjith  a  ratio 
t/a  =  1.02  (see  Figure  8),  which  is 
used  frequently  in  the  applications. 
The  following  figures  sho\;  the  radial 
and  tangential  stresses  Oj.  and  (Jq 
in  the  joint  for  different  materials. 
The  stresses  Qj.  are  all  compression 
and  ag  are  all  tension.  The  shear 
stress  Tj.0  can  be  considered 
zero  for  practical  purposes.  Thus, 
the  stress  trajectories  of  tension 
are  approximately  concentric  circles, 
and  the  compression  trajectories  are 
radial  straight  lines.  The  stress 
diagrams  shown  in  the  figures  corres¬ 
pond  to  the  following  cases: 


Figure  8.  Dimensions  and 
Coordinate  Axis. 


Figure 

No, 


Material 


9 

Isotropic  (steel) 

10 

Unidirectional 

glass 

composite  with 

Vp 

11 

Unidirectional 

glass 

composite  with 

Vp 

12 

Unidirectional 

glass 

composite  with 

Vp 

13 

Unidirectional 

boron 

composite  with 

Vp 

0.6 


0.7 


0.8 


0.7 


The  term  unidirectional  composite  indicates  that  all  the  fibers  in  the 
part  being  studied  (Figure  8)  are  parallel  to  the  circular  boundaries. 
Vp  indicates  the  volume  content  of  fibers. 


The  stresses  plotted  in  Figures  9  through  13  correspond  to  a  total  applied 
tension  force  P  =  2t  ;  in  other  words,  the  average  stress  (00)9=0 

unity.  From  the  analysis  of  these  stress  diagrams,  it  is  possible  to 
draw  the  following  conclusions: 


1.  The  peaks  of  both  stresses  and  Og  are  on  the 

internal  boundary. 

2.  The  maximum  difference  between  the  peak  values  of 
Cj.  is  less  than  157.  (from  1.14  to  1.32). 

3.  The  difference  in  the  stresses  Og  on  the  boundary 
in  the  critical  area  20‘’<6<80°  is  less  than  107.. 

4.  The  shapes  of  the  stress  diagrams  are  essentially 
the  same . 
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I 


I 


I 


Figure  9.  Stress  Diagrams  for  Isotropic  Material  (Steel). 
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a)  Radial  stress 


Load 

Direction 


0 


b)  Tangential  stress  oq 


Figure  11.  Stress  Diagrams  for  Glass  Composite  Material  (V 


Load 

Direction 


a)  Radial  stress  a^. 


b)  Tangential  stress  a, 

H 


gure  13.  Stress  Diagrams  for  Boron  Composite  Material  (V 


IL  is  therefore  possihli'  to  analyze  the  stresses  in  this  type  of  joint  by 
assuming  the  isotropy  of  the  material.  The  lesiiltant  stresses  will,  for 
design  juirposes,  be  sufficiently  accurate. 

Anotliei  type  of  parametric  study  was  performed  in  order  to  establish  the 
variation  of  the  stresses  with  respect  to  the  changes  of  the  ratio  t/z  . 
Two  types  of  materials,  isotropic  and  glass  composite  (Vp  =  0.7)  ,  were 
also  considered.  Table  1  gives  tlie  characteristics  of  the  joints  for 
which  stresses  are  shown  in  Figures  14  through  19. 


TABLE  I.  CHARACTERISTICS  OF  JOINTS 

Figure 

No. 

Material 

Ratio  — 
a 

Maximum, 

‘^r 

^  0  where 

Oj.  is  maximum 

14 

Isotropic 

0.57 

0.68 

1.07 

15 

Glass  composite 

0.57 

0.59 

0.98 

16 

Isotropic 

1.14 

1.28 

1.08 

17 

Glass  composite 

1.14 

1.22 

1.02 

18 

Isotropic 

1.71 

2.00 

1.29 

19 

Glass  composite 

1.71 

1.92 

1.20 

The  results  shown  in  this  table  confirm  the  conclusion  of  the  preceding 
parametric  study  referring  to  the  small  influence  of  the  material  proper¬ 
ties  on  the  maximum  stresses. 


Load 

Direction 


a)  Radial  stress 


Figure  14,  Stress  Diagrams  for  Isotropic  Material  (Steel). 


2 


Figure  15.  Stress  Diagrams  for  Glass  Composite  Material  (V 
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Direction 


Figure  16.  Stress  Diagrams  for  Isotropic  Material  (Steel) 


Load 

Direction 


a)  Radial  stress 


Figure  17.  Stress  Diagrams  for  Glass  Composite  Material 
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Figure  19.  Stress  Diagrams  for  Glass  Composite  Material 


lu'  i  iig 


'^r  (max)/'^avg  -’K'*  t/a  ,  -  ny^ 

till'  total  ajipl  it’ll  tonsilf  load. 

This  relation  can  be  expressed  with  minimal  error  by  the  linear  formula 

^r(max)  “  1.13  “  (54) 

Figure  21  shows  the  curve  of  “^th  t/a  as  the  independent  vari¬ 

able.  The  value  cjq  is  corresponding  to  the  point  in  which  Oj-  is  a  maxi¬ 
mum.  This  point  is  where  the  factor  of  safety  is  minimum,  as  we  will  see 
later.  The  curve  of  Figure  21  can  be  approximately  expressed  by  the 
following : 


Figure  20  is  a  plot  of  the  ratio 
the  average  stress  P/2t  and  P 


[l  +  0.079  (i)’] 


avg 


(55) 


Thickness  t/a 

Figure  20.  Maximum  Radial  Stress. 


Figure  21.  Value  of  cJq  Where 


Oj.  Is  Maximum. 


STRUCTURAL  DESIGN  OF  BOLT  JOINTS 


Failure  Criteria 


Once  the  stress  distribution  is  known,  it  is  necessary  to  adopt  a  failure 
criterion  in  order  to  establish  the  factor  of  safety  of  the  joint.  A 
failure  theory  widely  used  is  the  so-called  interaction  formula® 


(56) 


30 


t 


where  a„  ,  a,,  ,  and  T  are  the  macrotnechanicai  stresses  and  c  , 

^  y  j  ^  X  ^ 

Gy"  ,  and  '^xy  failure  stresses.  The  normal  stresses  and 

a  are  positive  (negative)  corresponding  to  tens/ on  (compression).  The 
failure  stress  Gx*  must  be  the  failure  stress  in  tension  if  G^^  is 
positive,  and  the  failure  stress  in  compression  if  Gj^  is  negative.  The 
same  criterion  is  used  with  respect  to  Gy  and  Gy*  .  In  the  formula, 

Gjj*  and  Gy*  must  always  be  positive.  The  interaction  formula  as  a 
failure  criterion  was  introduced  in  Reference  7.  Other  failure  criteria 
are  presented  in  References  11  and  12. 

In  the  following,  we  will  adopt  the  interaction  formula  (56)  as  the  repre¬ 
sentative  failure  criterion.  By  referring  equation  (56)  to  the  polar 
coordinates  r  and  0  ,  and  by  taking  into  account  that  Tr9  “  0  for  the 
critical  area  of  the  joint,  we  get  the  factor  cf  safety  p.: 


If  the  joint  is  made  in  the  manner 
indicated  in  Figure  22,  the  stress 
Oj.*  is  the  tension  failure  stress 
for  the  unidirectional  composite 
in  the  fiber  direction,  and  Gq* 
is  the  compression  failure  stress 
in  the  normal  direction. 

Design  Procedure 

The  design  procedure  for  a  bolt 
joint  will  be  illustrated  on  a 
joint  selected  from  the  Mohawk 
plane  of  Grumman  Aircraft  Company. 
This  joint  corresponds  to  the 
wing-fuselage  connection,  and  the 
ultimate  load  in  tension  is 
P  =  30,000  lb.  Information  about 
the  joint  design  is  given  in 
References  13,  14,  and  15. 


Fibers 


Figure  22.  Unidirectional  Com¬ 
posite  Joint. 


Figure  23  shows  the  preliminary  adopted  dimensions  by  using  glass  fibers 
and  epoxy  resin  with  a  fiber  volume  content  of  70%  to  757,. 
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Figure  23,  Dimensions  of  the  BoLt  Joint  (inch). 
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The  ratios  t/a  and  a.„„ 

a 

°^avg 

Then,  from  equations  (54)  and  (55),  we  get 

(^r^max  *  1*89  ^avg  »  “  1.22  ^ayg 

or,  with  the  value  of  <^avg  obtained  before, 

(ar)max  ^  75,600  psi  ,  00  =  48,800  psi 

The  failure  stresses  for  the  ui\idirectional  composite  employed’  are 
Gj.*  =  30,000  psi  ,  00*  =  225,000  psi 


are 


1.000  -  0.375 
0.375 


-  1.67 


30.000 


0.6  •  0.625 


40,000  psi 


L’he  failure  stress  in  compression  corresponds  to  the  0.27.  limit  of  strain. 
It  is  not  an  actual  'failure  stress."  By  using  this  value  and  on  the 
basis  of  equation  (57),  a  factor  of  safety  of  [i  =  0.38  results  and  the 
failure  load  will  be 


P*  =  0.38  •  30,000  =  11,400  lb 


However,  the  failure  load  obtained 
from  tests  is  about  26,000  lb,  as 
shown  in  the  discussion  of  experi¬ 
mental  results  in  this  report. 

This  discrepancy  between  analytical 
and  experimental  results  stems  from 
the  following  fact:  When  the  load 
reaches  a  value  of  12,000  lb,  a 
plastic  area  appears  as  shown  in 
Figure  24,  and  the  load-elongation 
diagram  presents  "jumps"  which 


Figure  24.  Plastic  Region, 


/ 
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indicate  the  successive  progress  of  the  plasticized  area.  Thus,  the  load 
obtained  from  the  linear  stress  analysis  only  indicates  the  beginning  of 
the  formation  of  the  plastic  areas. 

By  taking  bushings  of  larger  external  diameter,  we  can  make  joints  of 
smaller  ratio  t/a  .  The  following  two  values  are  considered: 

-  =  1.0  ,  -  =  0,66 

a  a 

The  corresponding  values  of  obtained  by  following  the  same  procedure 
as  above  result  in  the  values  in  Table  II. 


TABLE  II.  SAFETY  /ACTORS 

t 

a 

1  .0 

0.66 

‘^avg 

50,000 

62,500 

^‘^r^max 

DO , 500 

46,600 

^9 

53,900 

64,700 

u 

0.504 

0.582 

This  shows  that  the  factor  of  safety  increases  when  t/a  decreases. 
However,  in  practical  applications  of  joints,  other  limitations  exist  that 
do  not  permit  the  use  of  a  very  low  value  of  t/a  ,  A  ratio  t/a  =  I 
appears  to  be  quite  adequate  for  most  of  the  bolt  joints  of  composite 
materials. 

The  ultimate  failure  load,  which  is  the  maximum  load  that  can  possibly  be 
applied  to  the  joint  before  it  breaks  completely,  is  approximately  double 
the  load  given  by  the  described  procedure.  This  assumption  proved  to  be 
valid  only  for  0.9  ^  t/a  <  1.1  . 

EXPERIMENTAL  RESULTS 


The  experimental  program  was  divided  into  two  parts,  one  concerning  a 
simple  loop  joint  and  the  other  the  bolt  joint. 


34 


Loop  Joint 

An  experimental  program  concerned  with  the  stress  distribution  and  failure 
Loads  of  loop  joints  was  developed. 

The  joints  tested  to  failure  had  the  dimensions  indicated  in  Figure  25. 


n 


>—0,35  in. 


I II . 


0.70  in. 


I  igurc  25.  Dimensions  of  the  Loop  Joint. 

The  spcH  inuMis  i  ic  madi'  of  glass  composite  material  and  were  tested  in  ten¬ 
sion.  Till'  1  o.K.  v.is  apiilied  by  using  a  fixture  schematically  shown  in 
Figure  26.  A  (iiishing  was  not  used  between  bolt  and  joint,  and  special  care 
was  laki-n  to  msme  perfect  initial  contact  between  bolt  and  joint  at  the 
surface  of  load  transference.  1  n 


Figure  26.  Attachment  Device  for  Testing 


The  failure  loads  obtained  from  the  tests  are  summarized  in  Table  III, 


TABLE  III,  FAILURE  LOADS 

Dimensions  of  a 

Failure 

Avg 

No 

Simple  Cross  Section, 

load , 

Failure  Stress, 

in. 

lb 

psi  X  10  ^ 

1 

0.342  X  0,308 

12,900 

61.5 

2 

0.346  X  0.295 

11,250 

55.2 

3 

0,344  X  0.282 

10,400 

53.7 

h 

0,348  X  0.301 

13,350 

63,2 

5 

0.346  X  0,378 

13,100 

52,0 

The  ultimate  strength  of  the  joint's  material  in  simple  tension  was 
150-180  kpsi.  Therefore,  the  efficiency  of  the  loop  joint  is  about  0.33. 
In  other  words,  the  ultimate  strength  of  the  loop  joint  is  one-third  of 
the  strength  of  the  same  material  in  simple  tension.  This  conclusion  is 
limited  only  to  joints  with  a  ratio  t/a  =  1, 

Figures  27,  28,  and  29  illustrate  the  locati 
and  form  of  failure.  The  regions  of  resin  u 
failure  (opaque  areas  in  the  photographs)  en 
pass  three-fourths  of  the  total  section.  It 
also  possible  to  observe  that  the  failure  is 
produced  for  20°  <  0  <  80°  . 

A  failure  test  was  performed  by  using  a  plas 
homogeneous,  brittle  material  loop  joint. 

Figure  30  shows  the  failure  produced. 

Photoelastic  studies  were  conducted  to  estab 
lish  the  stress  distribution  in  both  ortho¬ 
tropic  and  isotropic  loop  joints.  Figure  31 
shows  the  isoclines  for  0°,  30°,  and  60°  for 
the  isotropic  joint.  Similar  patterns  were 
obtained  on  orthotropic  loop  joints  by  using 
the  reflection  method  of  photoelasticity. 

The  position  of  the  inclines  confirms  the 
analytical  results  in  the  sense  that  the 
shear  stresses  T^g  are  very  small.  In 
fact,  the  isocline  lines  are  approximately 
radial  straight  lines,  forming  angles  of 
0° ,  30°,  and  60°,  respectively. 


36 


Figure  28,  Loop  Joint  of  Composite 
Material  After  Testing 
to  Failure. 


Figure  29.  Loop  Joint  of  Composite 
Material  After  Testing 
to  Failure, 
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Figure  30.  Failure  of  an  Isotropic  Joint  of 
Brittle  Material. 


(a)  Isoclines  of  0°  (b)  Isoclines  of  30° 


I  » 


(c)  Isoclines  of  60° 


Figure  31.  Isoclines  for  0°,  30°,  and  60° 
for  the  Isotropic  Joint. 
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Bolt  Joint 


The  tested  bolt  joints  were  similar  to  those  shown  in  Figure  32,  Table  IV 
indicates  the  characteristics  and  ultimate  load  of  the  tested  bolt  joints. 


(a)  Composite  Joints 


(b)  Composite  Joints 


Bolt  Joints  Similar  to 
Those  Tested, 


\ 


TABLE 

IV,  SPECIMEN  CHARACTERISTICS  AND  ULTIMATE  LOADS 

1  _ _ _ ] 

Specimen 

No. 

Ratio 

t/a 

Ultimate 
Load ,  lb 

Figure 

No. 

Ob servations 

1 

1.67 

33 

-- 

2 

1.67 

34 

-- 

3 

1.67 

24,500 

-- 

-- 

4 

1.67 

-- 

-- 

5 

1.67 

26,700 

-- 

-- 

6 

1.67 

25,900 

-- 

-- 

7 

1.00 

45,400 

35 

Laterally  Supported 

8 

1.00 

44,500 

-- 

Laterally  Supported 

9 

1.00 

39,400 

36 

-- 

10 

1.00 

38,800 

-- 

-- 

11 

0.66 

42,200 

37 

-- 

12 

0.66 

41,500 

-  “ 

Specimens  7  and  8  were  tested  with  a  device  that  restricted  lateral  expan¬ 
sion  of  the  joint  under  high  loads.  In  this  way,  an  increment  of  the 
ultimate  load  of  10%  was  obtained. 


Figure  33.  Composite  Joint  After  Failure, 


40 


Joint  After  Failure 


Figure  38  represents  a  typical  load -displacement  curve  obtained  during  the 
tension  test  of  the  bolt  joints.  It  can  be  observed  that  there  are  several 
"jumps"  which  correspond  to  several  steps  of  failure  that  occurred  during 
the  test.  The  ultimate  load  shown  in  Table  IV  is  the  highest,  after  which 
the  joint  breaks. 


Elongation,  7. 


Figure  38,  Testing  Machine  Load-Displacement 
Plotted. 
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Figure  39  shows  the  specimen 
mounted  in  the  testing 
machine  for  tension  testing. 

FABRICATION  PROCESS 

The  following  procedure  was 
adhered  to  in  the  fabrication 
of  test  specimens  for  this 
program. 

Material 


1.  12  End  Roving  5994 

Glass 

2.  Resin  System 

a.  Epon  828  50  phr 

b.  Epon  1031  50  phr 

c.  NMA  (methyl 
endomethylene 
tetrahydro  phtha- 
lic  anhydride 
catalyst)  90  phr 

d.  Benzyldimethyla- 
mine  (BDHA  accel¬ 
erator)  0.50  phr 

Tooling 


Figure  39.  Specimen  Mounted  in  the 
Test  Machine. 


Winding  Tool  (  igure  40),  Aluminum-Constructed  Material 


Equipment 

1.  Temperature  Controlled  Hot-Melt  Coater 

2.  Fiber  Tensioning  Device 

3.  Winding  Machine  (special  throat-deep;  see  Figure  41) 

a.  Traverse-Controlled 

b.  Revolution  Infinite  Control 

Winding  Procedure 

1.  Set  winding  traverse  to  0.600  in.  to  yield  9  revolutions 
per  layer. 

2.  Set  fiber  tensioning  device  at  2-pound  fiber  tension. 

3.  Regulate  hot-melt  coater  temperature  to  170°-180°F. 

4.  Preheat  winding  fixture  to  160-180°F. 


5.  Wind  568  revolutions  at  2-pound  fiber  tension. 

6.  Apply  first-stage  side  plates;  close  side  plates  to  within 
0.062  In.  to  stops. 

7.  Place  assembly  In  a  225°F  regulated  press  platen;  close  side 
plates  just  short  of  resin  gelation. 

8.  Raise  temperature  to  250°?,  reside  for  1  hour  at  250°F,  and 
cure  at  325°F  for  1  hour.  Allow  to  cool  under  pressure. 

9.  Disassemble  tool  and  remove  precured  first-stage  laminate 
from  tool. 

10.  Remove  excess  material;  sandblast  faying  surface  for  second- 
stage  winding. 

11.  Wind  500  revolutions  with  fiber  tension  reduced  to  1.5  pounds. 
Total  winding  revolutions  of  1035. 

12.  Repeat  operations  6  through  9  using  second-stage  side  plates. 

13.  Remove  part  from  tool  and  report  all  physical  properties. 

STRESS  CONCENTRATION  AROUND  CIRCULAR  HOLES 
IN  ORTHOTROPIC  PLATES 

General  Formulas 

The  technical  Importance  of  the  stress  concentration  around  holes  Is  well 
known.  The  concentration  factors  have  been  established  for  the  most  Impor 
tant  cases  that  appear  for  Isotropic  materials  In  technical  applications. 
For  example,  In  a  tensioned  plate  with  a  circular  hole,  a  diagram  of 
stresses  appears  as  shown  In  Figure  42.  It  Is  possible  to  see  that 

the  stress  concentration  factor  Is  3,  because  the  peak  of  stress  Is  three 
times  larger  than  that  of  the  applied  stress. 


Figure  42.  Stress  Concentration  in  Isotropic  Sheet. 
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For  the  correct  use  of  structural  composite  materials,  it  is  important  to 
know  the  influence  of  the  relative  volume  and  directions  of  the  fibers  on 
the  concentration  factors.  This  is  of  fundamental  importance  in  the  case 
of  a  structure  with  alternative  loads,  because  the  maximum  stress  concen¬ 
tration  points  are  the  origin  of  the  failure  lines. 

In  this  part  of  the  report,  we  will  develop  the  analytical  elements  and  a 
computer  program  that  will  permit  us  to  find  the  concentration  factors  in 
an  orthotropic  plate  with  a  circular  hole.  As  external  loads,  we  consider 
tension  and  pure  shear.  By  combining  the  results  of  these  two  cases,  it  is 
possible  to  solve  almost  all  states  of  stresses,  assuming  that  the  hole's 
diameter  is  small  with  respect  to  the  dimensions  on  the  middle  plane  of  the 
plate. 

The  theoretical  elements  of  the  complex  solution  of  the  plane  state  of 
stresses  are  taken  from  Reference  16. 


The  fundamental  differential  equation  of  an  orthotropic  plate  in  plane 
stress  is 


a** 

bii  — +  (2bia  H- 

dy 


(58) 


where  bn  ,  bja.  bag  •  and  baa  are  the  elements  of  the  matrix  that 
connects  stresses  and  strains;  i.e., 


It  must  be  noted  that  the  matrix  b^j  is  the  inverse  of  matrix  [C]  men¬ 
tioned  earlier  in  this  report,  if  this  matrix  also  belongs  to  an  orthotropic 
plate. 


It  is  possible  to  rewrite  this  equation  in  the  following  form: 


with  the  constants  Oj 


and  Oa  defined  by 


OiOa 


ai  +  Oa 


2bi a  +  baa 

baa 


(59) 


(60) 
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In  other  words,  ax  and  aa  are  the  routs  of  the  equation  In  a 

haa  a  ■  (2bia  "h  ^>33)  a  +  bxx  =  0  (61) 

It  is  possible  to  demonstrate  by  energy  considerations  that  the  roots  of 
equation  (61)  are  real  and  positive  If  the  coefficients  bxi  ,  bia«  ^2a> 
and  b33  correspond  to  any  orthotropic  material. 

Taking  the  complex  variable  z  =  x  +  iy  and  its  conjugate  z  =  x  -  iy  , 
we  put  the  stress  function  in  the  following  form: 


i  =  f  (z  +  Yi  z)  +  g  (z  +  Yaz) 


(62) 


Substituting  this  function  into  the  differential  equation  (58),  we  can 
establish  that  the  real  part  of  equation  (62)  is  a  solution  if  Vi 
Ya  are  given  by 


Yi 


./or  -  1 

/oT  +  1 


Ya 


t/cta  ■  1 
/Oe  +  1 


(6?.) 


Because  tti  and  Oa  are  real  and  positive,  Yi  ^nd  Ya  between  -I 

and  1. 


Remembering  that  the  stresses  are 


d®* 

dx® 


xy 


y  $ 

9x3  y 


(64) 


By  substitution  of  equation  (62)  ,  we  obtain 


Ox  =  -  (1  -  Yi)  i"  (z  +  Yi^)  -  (1  -  Ya)  g"  (z  +  Yaz) 
Gy  =  (1  +  Yx  )*  f "  (z  +  Yxz)  +  (1  +  Ya)*  g"  (z  +  Yaz) 
Txy  =  -i(l  -  f"  (z  +  YxZ)  -  i(l  -  y|)  g"  (z  +  Yaz) 


(65) 


Introducing  the  function  t  =  '|'(x,y)  by 


9x9y 


9x" 


3  '  O'!  eta  ^  3 


9y- 


(1  -  Oi)(l  -  Og) 


9x9y 


(66) 
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the  displacements  are  given  by 


From  equations 

5\|f 

dx 


hy 


=  (bia  -  baa )  +  bi  1 

uy  -  (bia  -  baa)  §7  +  ^aa  gj 

(62)  and  (66),  it  is  possible  to  obtain 

4iVi  _  4iYa 

“  -  ~~y^  ■  "-'  v'a 

^Va 

=  rr77  s’ 


(67) 


(68) 


With  equations  (62)  and  (68),  the  displacements  given  in  equa¬ 
tion  (67)  become 

=  (1  +  Vi)(bia  -  aabaa)  ^  (z  +  Yiz) 

+  (1  +  Ya)(bia  -  ajbaa)  g'  (z  +  Vaz) 

Uy  =  i(l  -  YjXbia  -  Oiibaa)  f  (z  +  Yx^)  + 

+  i(l  -  Ya)(bia  *  ^a^az)  8*  (z  +  Ya^)  (69) 

Stress  function  (62)  must  produce  single-valued  stresses,  which  tend  to 
zero  at  infinity.  Also,  the  displacements  must  be  single-valued.  These 
conditions  are  satisfied  by  taking  as  function  i  the  real  part  of 

i  =  f(z  +  Yiz)  +  {A(1  -  YaXbia  -  Oabaa) 

+  iB(l  +  Ya)(bia  "  Ctibaa)}  (z  +  Yi^)  ln(z  +  Yi^) 

+  g(z  +  Yiz)  -  [A(1  -  YiXbia  -  ttibaa) 

+  iB(l  +  Yi)(bia  -  ttabaa))  (z  +  Ya^)  ln(z  +  Yaz)  (70) 

where  A  and  B  are  real  constants. 
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The  stresses  corresponding  to  the  function  (70)  are  obtained  by  using 
equations  (64 )  : 

=  -(I  -  Vi)^  (f"(z  +  Yiz)  +  [A(l  -  Y8)(bia  -  aabga) 

+  iB(l  +  YaXbia  -  a^baa)]  - - — =) 

z  +  Yi  zj 

-  (1  -  Ya)®Jg"(z  +  Yaz)  -  lA(l  -  Yi)(bia  -  O-ibaa) 

+  iB(l  +  Yi)(bia  -  ttabaa)]  - 

z  +  YazJ 

=  (1  +  Yi)®  +  YiZ)  +  [A(l  ^  Ya)(bia  ’  ttabaa) 

+  iB(l  +  YaXbia  -  ttibaa)]  - 

z  +  Yiz  J 

+  (1  +  Ya)*{g"(z  +  Yaz)  -  [A(l  -  Yi)(bia  -  ttjbaa) 

+  iB(l  +  YiXbia  -  ttabaa)]  "  -  j 


xy 


-i(l  -  y!)  |f"(z  +  Yiz)  +  (A(l  -  YaXbia  -  ttabaa ) 

+  iB(l  +  Ya)(bia  -  CXibaa)]  — 7^ — zl 

z  +  Yi  z  J 

-  (1  -  Ya)  |g"(z  +  YaZ)  -  tA(l  -  Yi)(bia  ‘  Otibaa) 


+  iB(l  +  Yi)(bi3  -  Qabaa)] 


Now  we  Introduce  the  two  new  functions  of  the  complex  variable 


z  =  re 


le 


(71) 


(72) 


defined  by  the  Integral  forms 


=  -  f 

2n  J 


2^  ze'^®  + 


-18 

ze  -  a 


('’r)r.a  <‘6 
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«  '  J, 


2n  -i8  , 
zc  +  a 


From  these  equations,  It  Is  possible  to  see  that  the  functions  P  and 
Q  are  related  to  the  tractions  on  the  boundary  of  the  circular  hole. 

In  these  equations,  0  is  the  angular  coordinate  of  the  point  on  the 
boundary  r=a  in  which  the  stresses  ffj.  and  are  taken.  The  func¬ 

tions  P  and  Q  are  finite  at  infinity  and 


{Re[P]}, 


-(Or), 


Substituting  equations  (71)  into  the  transformation  formulas 
Cj.  =  OjjCos^  6  +  Oy  sin*  0  +  2Tjjy  sin  0  cos  0 
Tj.0  =  (a  -  a^)  sin  0  cos  0  +  (cos*  0  -  sin*  0) 


we  obtain  the  polar  components  of  the  stresses: 


(2  -  y^zf 


f  (z  +  yj)  +  [A(i  -  Ya)(bia  -  agbaa) 


(z  -  Yaz)*  ( 

- = -  j  g" 

22  ( 


+  iB(l  +  Ya)(^ia  "  flibga)!  " '  '  _ | 

2  +  Yi*J 

'(2  +  Yaz)  -  tA(l  -  YiXbia  -  ftibaa) 
+  iB(l  +  YiXbia  -  Oebaa))  - 

Z  +  Y8*J 


i(z*  -  Y*1? 


^  I  f ' 


(z  +  Yiz)  +  [A(l  -  Ya)(bia  -  Oabaa) 


+  iB(l  +  Y3)(bia  -  dibaa)] 


-^1- 

+  YazJ 


(continued) 


1 


i  (z  -  \|z®)  f  _ 

*  - — -  Sk"^z  +  Yzz)  -  [A(l  -  YiXbis  -  ttibsg) 

+  iB(l  +  Yi)(b  la  ■  Q’abaa)]  ;  (76) 

z  +  YazJ 

By  putting  these  stresses  into  equations  (7A),  functions  P  and  Q 
become 

p  -  ^  .  .(v^) .  > 


+  Ad 


^  /l  +  Yi  Ayi  \ 

Ya)(bia  -  Q-abaa )  I  - -  -  “"^a  1 

^  ^  z 

-  18(1  +  VaXbi,  -  a,b„) 


-  A(l  -  Yi)(bia  *  ttibaa) 


iB(l  +  YiXbia  -  ttabaa) 


n  +  Ya 


^Ya 


z  +  -*- 


/l  -  Ya  ^Y8  \ 

9  9 


+  (A  +  iB)((l  +  Yj)(l  -  YaXbia  -  aabja) 


(I  -  YiXl  +  YaXbia  -  di^aa))  7  (77) 


Now  we  define  the  two  new  variables  §  and  T|  by  the  following  equations: 


x[‘X‘ 


AY^a* 


x[‘X 


I  - 


AYiaa 


.8.  h 


(78) 
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For  a  particular  problem,  functions  P  and  Q  can  be  computed  by  using 
equations  (7A),  because  (aj.)j._g  and  (^r0)r=a  known  from  tlie  bound¬ 

ary  conditions.  By  integration  of  tlie  stresses  on  the  boundary  of  tlie 
hole,  we  find  the  resultant  forces 

X  =  2nA(l  -  Vi)(l  -  Va)(ai  -  aa)baa 

Y  =  2ttB(1  +  YjXl  +  YaXtti  -  aa)baa  (81) 

Equations  (73)  can  be  used  to  find  V(z)  and  W(z)  ,  and  by  solving 

equations  (81),  constants  A  and  B  can  be  established.  Then,  from 
equations  (79)  and  (80),  we  obtain  f"(§)  and  g"  (^)  •  The  stresses  are 
given  by  equations  (71),  where  we  now  have 

§  =  z  +  YjZ  ,  T]  =  z  +  VaZ  (82) 

Stress  Concentration  for  Pure  Tension 

The  tension  load  T^^  is  applied  at  the  x  direction,  as  shown  in 
Figure  43. 


*x 

Figure  43.  Tension  Load  Applied 


at  the  X  Direction. 

If  there  is  no  hole  in  the  plate,  the  stresses  are 

^x  “  ^x  •  “  0  ,  'xy  “  ® 

or,  in  polar  coordinates,  according  to  the  formulas  (75), 

T 

Qj.  =  (I  +  cos  29) 
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=  -r-  (:  -  cos  26) 


T^e  =  -  —  sin  26 


It  is  now  necessary  to  find  a  stress  system  which  tends  to  zero  at  the 
infinite  and  which  takes  the  values 


(I  -t-  cos  26) 


^  (I  -  cos  26) 


i  i 


C^re^s  =  ‘f  26 


when  r  =  a 


Thus,  from  equations  (73),  we  get 


With  the  variables  ^  and  T]  introduced  by  equations  (78),  we  obtain 
the  following  from  equations  (84): 


P(5)  -  Q(§)  =  P(T1)  -  Q(T1)  “  T 


Tx  <1  -  I 


P(§)  +  Q(§)  = 


Avja- 


P(T1)  +  Q(T1)  = 


(  f  ^Yia*\V 

V. 

r— II 

f  4Yaa« 

Ya  U+  I  --4- 


Because  the  stresses  in  equation  (85)  do  not  produce  a  resultant  force  on 
the  hole,  constants  A  and  B  are  zero,  as  equations  (81)  demonstrate. 
Thus,  from  equations  (79)  and  (80),  and  by  taking  into  account  equations 
(77)  and  (87),  we  get 


f"(5) 


g"  (Tl) 


Tx(^  +  Va) 
''•Vi  (Yi  -  Ya) 


Tx(l  Yi) 
^Ya(Yi  -  Ya) 


7 


(88) 


Using  these  equations  and  remembering  equations  (82),  then  the  following 
stresses  result  from  equations  (71): 


(1+  Yi)(l  +  Ya)Tx  ,  Tx 
"  Ay^Ya  A(>^  -  Va) 


(1  +  Ya)(l  -  Yi) 

YiXi 

(1  +  Yi)(l  -  yJ“ 


YaXa 


(1  +  Yi)(l  +  Ya)Tx 


AYi  Ya 


A(Yl-  Ya) 


[  (1  +  Ya)(l  + 

\  YiXi 


Yi)’ 


(1  H-  Yi  )(1  +  Ya)* 


YaXa 


iT. 


xy  A(Y:  -  Ya) 


(1  +  >^)(i  -  Yf)  a 

YiXi 


+  Yi)(^  -  V^)  ] 

YaXa  J 


(89) 


where  the  following  notation  was  used: 

if, 

Xj  -  Rl^  e  , 


if  a 


Xa  =  Ra  e 


(90) 


with 


n 

Rj  e  1 


1  - 


Aa*Y, 


(z  +  Yi  z)* 


Ra  e 


ifa 


Aa®Ya 


(z  +  Yaz)* 


(91) 
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The  quantities  Rj  and  Rj  are  chosen  as  positive,  and  and 

are  to  lie  between  -n  and  tt  .  Also,  it  must  be  emphasized  that  the 

actual  stresses  are  the  real  part  of  complex  expressions  (89). 

Stress  Concentration  for  Pure;  Shear 

A  shear  load  Tj^y  is  applied  on  the  plate  with  the  circular  hole,  as 

shown  in  Figure  44.  If  the  hole  is  absent,  the  stresses  are 

=  0  ,  Oy  =  0  ,  Txy  =  Txy  (92) 

or,  in  polar  coordinates, 

^r  "  ^xy 

00  =  -Txy  sin  20 

"^re  =  cos  20  (93) 

xy 

X 


Figure  44.  Shear  Load  Applied  on 

Plate  With  Circular  Hole. 

Thus,  to  keep  the  hole  boundary  free  of  stresses,  it  is  necessary  to  sum 
to  the  stresses  (93),  the  stresses  which  take  the  values 

(^i-)s  “  "^xy  »  ("'^rQ^s  “  ""^xy  cos  20  (94) 

on  the  hole  and  zero  at  the  infinite. 
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By  urfin^  equations  (73),  we  obtain 


P  =  Q 


(95) 


and,  from  equations  (8-'t), 


P(!)  -  Q(!) 


P(?)  +  Q(§)  = 


P(11)  +  Q(T1)  = 


(96) 


(97) 


The  constants  A  and  B  of  equations  (81)  are  zero  because  the  resultant 
force  (X,  Y)  on  the  hole  is  zero.  With  equations  (96)  and  (97),  from 
equations  (79)  and  (80),  we  have 


f"(?) 

g"  (H) 


i  T 


EL 


2Vi  (Yi  -  Ya  ) 


( 


^Ya  ('fi  ■  Ya  ) 


(‘ 

(‘ 


4Yia®N^ 


(98) 


With  these  results  and  with  equations  (82),  from  equations  (71)  the 
stresses  are  tlie  real  part  of  tlie  expressions 


..  _  i 

(  (1  -  '^i)® 

(1  -  y-J 

2(Y:  -  V) 

1 

YaX, 

^  Ixy 

f(l  +  >i)" 

(1  +  Ya) 

"  ■  2(Vi  -  V,- 

\  iXj 

VaX',! 

(99) 
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(100) 


i 


xy 

The  meaning  of  Xj 
Numerical  Results 


ViVa  -  ^ 
^Vi  Va 

and  Xg 


is  given  in  equations  (90)  and  (91). 


By  using  equations  (89)  and  (99),  it  is  possible  to  study  the  stress  con¬ 
centration  around  a  circular  hole  on  an  orthotrojjic  plate  under  axial  and 
shear  loads,  A  computer  program  based  on  these  complex  expressions  was 
developed . 

A  representative  value  of  the  stress  concentration  originated  by  a  circular 
hole  is  the  stress  (Ox)ni  point  A  when  the  applied  load  produces 

unit  stress  in  the  plate  without  a  hole.  By  using  the  computer  program, 
we  obtained  the  results  plotted  in  Figure  45.  It  can  be  observed  that  the 
maximum  stress  concentration  factor  is  obtained  when  the  axial  load  is 
applied  in  the  direction  of  the  fibers  (x  ,  in  this  case). 


Figure  45,  Stress  Concentrations  (maximum  stresses  c^ra) 
in  a  Unidirectional  Composite  With  a  Circular 
Ho  1  e . 
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The  magnitude  of  the  concentration  factor  indicates  tlint  tliis  plienomenon 
is  more  critical  than  in  the  isotropic  inati>r  i  .a  I  s ,  where  the  factor  is 

only  3. 


It  must  be  noted  that 
material  constants  do 
the  hole  must  be  made 


this  analysis  is  based  on  tlie  liypothesis  that  the 
not  depend  on  x  and  y.  To  satisfy  this  condition, 
by  cutting  the  fibers  as  indicated  in  Figure  46. 


Figure  46.  Orthotropic  Sheet. 


A  more  rational  form  for  the  fibers  is  shown  in  Figure  47,  where  fibers 
are  not  cut.  The  corresponding  stress  analysis  cannot  be  performed  with 
expressions  (89)  and  (99)  because  now  the  elastic  constants  change  from 
point  to  point.  In  other  words,  the  material  is  not  quasi -homogeneous. 
However,  the  stress  analysis  of  this  case  can  be  made  by  using  the  method 
developed  on  page  46. 


Figure  47.  Nonorthotroplc  Sheet. 
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The  numerical  study  was  performed  for  unidirectional  composites,  in  which 
the  importance  of  the  stress  concentration  factor  is  more  evident. 
However,  formulas  (89)  and  (99)  can  also  be  applied  for  multilayered 
composite  materials  by  considering  the  adequate  constants  b^j  . 


CONCLUSIONS  AND  RKCOMMKNDATIONS 


The  introduction  of  localized  high  loads  into  a  thin  sheet  c.in  be  accom¬ 
plished  effectively  by  the  use  of  bolt  joints;  contrary  to  the  metal¬ 
lic  joints,  there  are  no  discontinuities  present  when  joints  made  of 
filamentary  composites  are  used.  A  furtliei  advantage  of  filamentary  joints 
is  their  behavior  when  loaded  to  failure.  For  example,  after  the  first 
break  occurs,  tlie  joint  can  still  be  loaded  to  almost  twice  the  failure 
load.  Since  the  first  break  is  very  easily  detectable,  it  is  possible  to 
avoid  continuation  of  a  catastrophic  failure  and  still  maintain  the  load- 
carrying  capability  of  a  structure. 

The  stress  analysis  of  composite  joints,  as  presented  herein,  shows  that 
an  effective  safety  factor  can  be  obtained  that  is  within  the  accuracy 
demanded  by  technical  requirements. 

In  a  circular  form  cutout  in  a  unidirectional  composite,  the  stresses  can 
reach  values  twice  as  high  as  thos?  attained  in  an  isotropic  material  if 
the  fibers  terminate  in  the  rim  of  the  cutout.  For  this  reason,  it  is 
recommended  that  the  fibers  be  led  in  a  continuous  manner  around  the  cutout 
as  far  as  fabrication  requirements  permit.  Impelled  by  the  great  poten¬ 
tials  of  composite  joints,  further  studies  of  other  shapes  of  joints  and 
types  of  loadings  are  reconmended. 
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